#市町村情報
shichoson_info <- st_drop_geometry(shichoson_geo)
gw_2_data <- st_transform(gw_2_data, st_crs(shichoson_geo))
intersections <- st_intersects(shichoson_geo, gw_2_data, sparse = FALSE) %>% which(. == TRUE)
shichoson_geo["passed"] <- 0
shichoson_geo[intersections,"passed"] <- 1
passed_pre_geo <- shichoson_geo[which(shichoson_geo$N03_001 %in% unique(shichoson_geo[intersections,1]$N03_001)),]
mapview(
passed_pre_geo,
zcol = "passed",
layer.name = NA,
col.regions = c("white", "red"),
popup = NA,
label = passed_pre_geo$N03_004,
legend = F
) +
mapview(gw_2_data,
label = NA)
track_points <- track_points %>%
mutate(lon = st_coordinates(.)[,1],
lat = st_coordinates(.)[,2])
track_points <- track_points %>%
mutate(dist = geosphere::distGeo(cbind(lon, lat),
cbind(lag(lon), lag(lat)))) %>%
mutate(dist = ifelse(is.na(dist), 0, dist)) %>%
mutate(cumdist = cumsum(dist))
ggplot(track_points, aes(x = cumdist / 1000, y = ele)) +
geom_line(color = "steelblue") +
labs(x = "Distance (km)", y = "Elevation (m)", title = "Elevation Profile") +
theme_minimal()
